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Abstract. We report 250 GHz (1.2 mm) observations of redshift 1.8 < z < 2.8 optically luminous (Mb < —27.0), 
radio quiet quasars using the Max-Planck Millimeter Bolometer (MAMBO) array at the IRAM SO-metre telescope. 
Nine quasars were detected and for 26 quasars Scr flux density limits in the range 1.8 to 4 mjy were obtained. 
Adopting a typical dust temperature of 45 K, the millimeter emission implies far-infrared luminosities of order 
10“ Lo and dust masses of ~ 10® Mq. Applying a statistical survival analysis to our total sample of 4S detected 
and 95 undetected quasars at 2 ; « 2 and 2 ^ 4, we find that there is no apparent difference in the far-infrared (FIR) 
luminosities, hence the star formation rates, of QSOs at 2 « 2 and at S.6 < 2 < 5. This differs from radio galaxies, 
for which the FIR luminosity was found to increase with redshift (Archibald et al. 2001). We furthermore find 
that there is no strong correlation between the far-infrared and optical luminosities, confirming previous results 
obtained on smaller samples. 

Key words. Galaxies: formation - Galaxies: starburst - Galaxies: high-redshift - Quasars: general - Cosmology: 
observations - Submillimeter 


1. Introduction 

The relation between the growth of the central black hole 
and the formation of the bulge stars is a key issue for 
the formation and evolution of galaxies. Probing star- 
burst activity in the host galaxies of high redshift quasars 
has therefore become a key area of observational cos¬ 
mology. Since the formation of massive stars often oc¬ 
curs in heavily obscured regions, star formation is best 
traced through the far-infrared (FIR) dust emission (e.g., 
Sanders & Mirabel 1996), the peak of which is red-shifted 
into the (sub)millimeter atmospheric windows for sources 
at redshifts 2 > 1. Deep (sub)millimeter blank field sur¬ 
veys using SCUBA (e.g., Scott et al. 2002) and MAMBO 
(Bertoldi et al. 2000a,b) have now detected over one hun¬ 
dred sources. For most of these objects optical obscura¬ 
tion prevents a spectroscopic redshift determination and 
any detailed study of the relation between star formation 
and Active Galactic Nuclei (AGN). Therefore, such stud- 
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ies must currently rely on pointed (sub)millimeter obser¬ 
vations of radio galaxies and of optically selected quasars. 

Observations at 1.2 mm and 0.85 mm using MAMBO 
and SGUBA, respectively, have lead to the detection of 
more than 60 high-redshift quasars (Omont et al. 1996 
and references therein for earlier work; Omont et al. 2001; 
Carilli et al. 2001a; Isaak et al. 2002; and present paper; 
Priddey et al. 2002), radio galaxies (Archibald et al. 2001), 
and X-ray selected AGN (Barger et al. 2001; Page et al. 
2001). Whether most of the dust FIR emission is due to 
heating by the UV/X radiation of the AGN or by mas¬ 
sive stars remains unclear. Only for a few cases there are 
searches for CO line emission deep enough for achieving 
detection and providing good evidence that a substantial 
fraction of the FIR emission must be caused by star for¬ 
mation (Guilloteau et al. 1999; Cox et al. 2002a,b; Carilli 
et al. 2002). 

Up to now most of the studies of the (sub)millimeter 
emission from high redhift QSOs have focused on sources 
at redshifts 2^4, for which 40 detections were reported 
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(Omont et al. 1996, Omont et al. 2001, Carilli et al. 2001a, 
Isaak et al. 2002, and references therein). In order to study 
the quasar population at the peak of their space density, 
which is at a later epoch than that probed in previous 
studies, we extended the millimeter observations of op¬ 
tically luminous quasars to z « 2. This paper reports on 
1.2 mm continuum observations of a sample of redshift 1.8 
to 2.8 QSOs with optical luminosities in excess of 10^^ Lq 
(Mb < -27.0). 

We adopt a A-cosmology with Hq = 65 km s“^ Mpc“^, 
Ha = 0.7 and flm = 0.3. However, for the definition 
of the rest-frame absolute B-band magnitudes, Mb, we 
have used the standard Einstein-de Sitter cosmology with 
Hq = 50 kms“^ Mpc“^ and qo=0.5, in order to ease the 
comparisons with standard quasar catalogues and lumi¬ 
nosity functions. Throughout this paper, the far-infrared 
luminosity (Tfir) is defined as the luminosity of a mod¬ 
ified black-body with a dust temperature and emissivity 


index as given in Sect. 3.1 


2. Source Selection and Observations 


Thirty-five optically luminous (Mb < —27.0) radio-quiet 
QSOs at redshifts 1.8 < z < 2.8 were observed from 60 
sources we had selected from the catalogue of Veron-Cetty 
& Veron (2000). As the result of our random selection 
among the sources listed by Veron-Cetty & Veron, the ob¬ 
served sources derive from a broad variety of QSO surveys. 
However, about half of them were originally identihed in 
the Hamburg Quasar Survey (Engels et al. 1998; Hagen; 
Engels & Reimers 1999, littp://www.hs.uni-hamburg.deD . 
The observed sources have a distribution in rest-frame ab¬ 
solute B-band magnitudes, Mb, between —27.0 and —29.5, 
which is comparable to the Mb distribution of z ^ 4 PSS 
quasars studied by Omont et al. (2001) - see Sect. For 
the present sample at z 2, the values of Mb were deter¬ 
mined including J and H magnitudes obtained from the 
2MASS near-infrared survey - see Priddey et al. (2002). 
However, the Mb values are more uncertain for the z ^4 
sample because 2MASS data are lacking for most sources. 
The selected QSOs are at high declination (> 15° and 
mostly > 35°) to favor observations done at low air mass. 
The list of the observed sources is given in Tables 1 and 
2 . 


The observations were made during the winter of 
2000-2001 using the 37-channel Max-Planck Millimeter 
Bolometer (MAMBO; Kreysa et al. 1999a, 1999b) array 
at the 30-meter IRAM telescope on Pico Veleta (Spain). 
The sources were observed with the array’s central chan¬ 
nel, using the standard on-off mode with the telescope 
secondary chopping in azimuth by 50'' at a rate of 2 Hz. 
The target was positioned on the central bolometer of 
the array, and after 10 seconds of integration, the tele¬ 
scope was nodded so that the previous off beam became 
the on beam. A typical scan lasts for 12 or 16 such 10 
seconds sub-scans. The pointing was checked frequently 
on nearby continuum sources, and was found to be sta¬ 
ble within ^ 2". The sky opacity was measured regularly, 
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Fig. 1. Observed 1.2 mm flux density versus Mb, the op¬ 
tical absolute magnitude in the rest-frame B band, of the 
z Ks 2 QSOs. The sources detected at 1.2 mm (Tablel) 
are shown as filled symbols. The open symbols are upper 
limits for sources not detected at 1.2 mm (Table2): for 
sources with negative signal, a 1 u upper limit is plotted; 
for sources with positive signal, signal -I- Ict is plotted as 
an upper limit. 


with zenith opacities at 1.2 mm varying between 0.08 and 
0.4. Gain calibration was performed using observations of 
Mars. We adopted a calibration factor of 12500 counts per 
Jansky, which we estimate to be reliable to within 20%. 
The point-source sensitivity of MAMBO during the ob¬ 
servations was ~ 30mJys^/^. 

The data were analyzed using the MOPSI software 
(Zylka 1998). Correlated sky-noise was subtracted from 
each channel; it was computed for each channel as 
a weighted mean of the signals from the eight best- 
correlating surrounding channels. Skynoise subtraction re¬ 
duces the noise in the signals by typical factors of 2 to 3. 
For each quasar, the total on-target plus off-target observ¬ 
ing time was always greater than 650 sec. The mean r.m.s 
noise of the coadded signals is « 0.8 mJy. As for Omont et 
al. (2001), we defined a quality factor for the detections (A 
for good and B for poor) based on the consistency of the 
flux densities measured on different dates, on the stability 
of the pointing, on calibration uncertainties and general 
observing conditions (see Table 1). 


3. Results 

3.1. General properties 

Thirty-five z ~ 2 quasars were observed at 1.2 mm. Nine 
quasars were detected at levels > 3 cr, of which seven have 
signals ^ 4cr, where a w 0.8 mJy. The 1.2 mm flux den¬ 
sities of these quasars are listed in Table 1. Twenty-six 
sources were not detected with 3 a flux density upper lim¬ 
its in the range 1.8-4 mJy (Table 2). The weighted mean 
signal of the non-detections is 0.71 ± 0.16 mJy. 
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Table 1. Quasars detected at 1.2 mm. 


Source 

z 

Mb 

R.A. Dec. 

(J2000.0) 

Flux Density 
[mJy], ± 1(7 

time 

[sec] 

Quality 

1.4 GHz 
[mJy], ± Ict 

KUV 08086+4037 

1.78 

-27.0 

08 12 00.5 

+40 28 14.0 

4.3 ±0.8 

1336 

A 

- 

[VV96] J093750.9+730206 

2.52 

-28.5 

09 37 48.7 

+73 01 58.0 

3.8 ±0.9 

665 

B 

* 

HS 1002+4400 

2.08 

-28.3 

10 05 17.5 

+43 46 09.0 

4.2 ±0.8 

1654 

A 

- 

HS 1049+4033 

2.15 

-28.2 

10 51 58.6 

+40 17 36.0 

3.2 ±0.7 

1818 

A 

- 

[VV96] J110610.8+640008 

2.19 

-29.3 

11 06 10.8 

+64 00 08.0 

3.9 ± 1.1 

2228 

A 

* 

[VV96] J140955.5+562827 

2.56 

-28.4 

14 09 55.5 

+56 28 27.0 

10.7 ±0.6 

1832 

A 

- 

[VV96] J154359.3+535903 

2.37 

-28.3 

15 43 59.3 

+53 59 03.0 

3.8 ± 1.1 

992 

B 

- 

HS 1611+4719 

2.35 

-27.7 

16 12 39.9 

+47 11 58.0 

4.6 ±0.7 

3319 

A 

- 

[VV2000] J164914.9+530316 

2.26 

-28.2 

16 49 14.9 

+53 03 16.0 

4.6 ±0.8 

1655 

A 

0.70 ± 0.09 


NOTE - The radio flux densities at 1.4 GHz are extracted from the VLA FIRST survey. The symbol * indicates that 
there is no data available, and the minus symbol indicates a 6 a upper limit to the flux density of 1 mjy. 


Fig.g displays the observed 1.2 mm flux density S '250 
versus Mb, the optical absolute magnitude in the B rest- 
frame band, which represents the bolometric luminosity, 
Lboi- Leaving aside J1409, which has an exceptionally high 
1.2 mm flux density of 10.7 ±0.6mJy, the other eight de¬ 
tected quasars have comparable flux densities of 3-4 mJy. 
There is no obvious correlation between Mb and 5'25o, a 
point which will be discussed in Sect. 

With a single measurement at 250 GHz (at z ~ 2.2, 
this corresponds to an emitted frequency of ^ 800 GHz or 
a wavelength of 375/rm), it is not possible to distinguish 
whether the observed millimeter emission is thermal dust 
or synchrotron emission. The VLA 1.4 GHz FIRST survey 
(Becker, White, & Helfand 1995) provides a first estimate 
on the radio to millimeter slope. No radio source within 
30" of six of the z w 2 detected quasars is found to a 
6 a limiting flux density of w 1 mJy. Two mm-detected 
QSOs are not covered by the FIRST survey, and J1649 
coincides with a radio source with a 1.4 GHz flux density 
of 0.70 ± 0.09 mJy (Table 1). For the quasars not detected 
at 1.2 mm (Table 2), two sources have a radio counter¬ 
part with 1.4 GHz flux densities of 2mJy. In all cases, 
the 1.4 GHz flux densities are too small to consider the 
quasars as radio loud. For the detected source J1649, the 
millimeter-radio spectral index is « 0.37, well above the 
extrapolation of a synchrotron spectrum and close to the 
value expected at z ^ 2 for starbursts (e.g., Carilli & Yun 
1999). Its case is similar to the QSOs PSS J1048-I-4407 
and PSS J1057-I-4555 detected at 1.2 mm and 1.4 GHz 
with comparable flux densities (Omont et al. 2001; Carilli 
et al. 2001b). We shall assume in the following that the 
millimeter flux density of the z « 2 quasars detected at 
1.2 mm is thermal in nature, as supported by their very 
recent deep VLA 1.4 GHz observations (Petrie et al. in 
preparation). 

For consistency with our earlier treatments (Omont et 
al. 2001), we adopt a dust emissivity index of /3 = 1.5, 
and a dust temperatures of 45 kU In the redshift range 

^ Adopting different values for the dust temperature and /3, 
e.g., 40 K and 2.0 (Priddey & McMahon 2001), will not change 
the main conclusions of this paper 


2 < z < 3, the FIR luminosity can then be expressed as a 
function of the 1.2 mm (250 GHz) flux density, S'250, as 

Lfir ~ 4.7 X IO 12 (S 25 o/niJy) Lq . (1) 

The corresponding proportionality factor for the z ^ 4 
quasars is 3.5 x 10^^. The derived FIR luminosities for the 
z « 2 quasars detected at 1.2 mm are « 10^^ L©, with 
corresponding dust masses of a few 10® x (7.5/Kd23o) M©, 
where Kd 230 is the dust absorption coefficient at 230 /im 
in units of cm^g“^. 

Like for the QSOs detected at z ^ 4, for the newly 
detected objects Lfir is about one tenth of the optical 
luminosity Lopt- If a substantial fraction of Lfir arises 
from young stars, such luminosities imply star formation 
rates approaching 10® M© yr“®, comparable to the values 
derived for z ^ 4 quasars (Omont et al. 2001; Carilli et al. 
2001a). 

3.2. Comments on Individual Sources 

- VV96 J140955.5-b562827 (SBS 1408-b567) 

This quasar has an exceptionnally strong 1.2 mm flux 
density which puts it among the seven known high 
z sources with S 250 ^ 10 mJy. Four of the latter are 
known to be strongly lensed, while there is presently no 
indication of a strong amplification for SBS 1408-1-567. 
The optical spectrum of SBS 1408-1-567 shows spec¬ 
tacular detached broad absorption lines (BAL) with 
double troughs (Korista et al. 1993). Deriving the red- 
shift from the unpublished Keck HIRES spectrum of 
Barlow & Junkkarinen (1994) - and Junkkarinen (pri¬ 
vate communication), we recently detected CO in this 
QSO using the IRAM Plateau de Bure interferometer 
(Beelen et al. in preparation). 

- VV96 J154359.3-b535903 (SBS 1542-f541) 

This source was discovered in the Second Byurakan 
Survey (Stepanyan et al. 1992) and has many interest¬ 
ing properties: BAL with a very high ionization degree 
(Teller et al. 1998), associated absorption system and 
damped Lya (DLA) absorption system, and a strong 
X-ray absorption (Green et al. 2001). 
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Table 2. Quasars with upper limits at 1.2 mm. 


Source 

z 

Mb 


R.A. Dec 

(J2000.0) 


Flux Density 
[mjy], ± 1(7 

time 

[sec] 

1.4 GHz 
[mJy], ± 1(7 

HS 0749+4259 

1.90 

-28.9 

07 

50 

54.7 

+42 

52 

19.0 

1.3 ±0.9 

1652 

- 

HS 0800+303H 

2.02 

-28.6 

08 

03 

42.1 

+30 

22 

54.0 

1.8 + 1.1 

824 

- 

HS 0808+1218''‘ 

2.26 

-27.8 

08 

10 

57.0 

+12 

09 

14.0 

1.1 ±0.9 

1583 

- 

HS 0830+18331' 

2.27 

-28.1 

08 

32 

55.7 

+18 

23 

01.0 

1.3 ±0.9 

1414 

- 

HS 0834+1509t 

2.51 

-28.1 

08 

37 

12.6 

+14 

59 

17.0 

0.6 ±0.6 

2391 

- 

SBSS 0910+586 

1.95 

-27.1 

09 

14 

25.8 

+58 

25 

19.0 

0.4 ±0.6 

839 

* 

[VV96] J092230.1+710937 

2.43 

-27.1 

09 

22 

30.1 

+71 

09 

37.0 

1.5 ±0.8 

1666 

* 

HS 0932+2410 

2.30 

-27.1 

09 

35 

34.0 

+23 

57 

20.0 

1.4 ±0.6 

2074 

- 

[VV96] J093935.1+361001 

2.03 

-27.1 

09 

39 

35.1 

+36 

40 

01.0 

1.5 ± 1.0 

1498 

- 

[VV96] J095845.5+470324 

2.48 

-27.7 

09 

58 

45.5 

+47 

03 

24.0 

0.7 ±0.7 

1647 

- 

HS 1110+3813 

2.29 

-28.1 

11 

12 

51.0 

+37 

57 

15.0 

0.0 ± 1.3 

833 

- 

KUV 11467+3725 

2.22 

-27.8 

11 

49 

20.2 

+37 

08 

29.0 

1.0 ±0.9 

1667 

- 

[VV96] J121010.2+393936 

2.40 

-27.7 

12 

10 

10.2 

+39 

39 

36.0 

0.9 ±0.8 

1492 

- 

[VV96] J121303.1+1714221' 

2.54 

-28.0 

12 

13 

03.1 

+17 

14 

22.0 

-1.9 ± 1.5 

656 

1.97 + 0.16 

[BCF91] 524 

2.85 

-28.0 

13 

04 

12.0 

+29 

53 

49.0 

0.2 ± 1.0 

329 

- 

[BBE90] 130623+283002 

2.21 

-27.9 

13 

09 

17.2 

+28 

14 

04.0 

-1.6 ± 1.0 

1648 

- 

[VV96] J140148.4+543859 

2.37 

-27.3 

14 

01 

48.4 

+54 

38 

59.0 

0.3 ±0.9 

842 

- 

SBSS 1417+596 

2.31 

-27.6 

14 

19 

06.4 

+59 

23 

12.0 

-3.1 ± 1.6 

831 

* 

[VV96] J160637.6+173516 

2.32 

-27.4 

16 

06 

37.6 

+17 

35 

16.0 

1.6 ±0.7 

2176 

- 

HS 1616+3708 

2.49 

-27.4 

16 

18 

15.5 

+37 

01 

03.0 

1.0 ±0.6 

2831 

- 

[VV86] J162645.7+642654 

2.32 

-28.8 

16 

26 

45.7 

+64 

26 

54.0 

0.7± 1.1 

3470 

* 

HS 1707+4602 

2.29 

-27.6 

17 

09 

04.9 

+45 

59 

08.0 

2.1 ± 1.5 

2487 

- 

[VV96] J171635.4+532815 

1.94 

-28.8 

17 

16 

35.4 

+53 

28 

15.0 

1.3 ± 1.1 

2469 

1.90 + 0.16 

HS 1754+38181' 

2.16 

-27.9 

17 

56 

39.6 

+38 

17 

52.0 

1.1 ± 1.1 

3282 

* 

[VV96] J183825.0+510558 

1.98 

-29.3 

18 

37 

25.3 

+51 

05 

59.0 

-1.0 ±0.8 

2452 

* 

HS 2134+15311' 

2.13 

-28.1 

21 

36 

23.7 

+15 

45 

08.0 

-2.7 ± 1.4 

2494 

* 


NOTE - The radio flux densities at 1.4 GHz are extracted from the VLA FIRST survey. Symbols as for Table 1. 
^ Quasars observed (and not detected) at 850 fim by Priddey et al. (in preparation). 


- VV96 J110610.8+640008 (HS 1103+6416) 

It is the brightest optical source of the sample. Its 
spectrum was analysed by Koehler et al. (1999), us¬ 
ing a high resolution HIRES/Keck optical spectrum 
and a UV HST spectrum. The most striking feature is 
a strong, complex absorption system at z = 1.9. 

- KUV 08086+4037 

The only reference in SIMBAD is the discovery pa¬ 
per (Darling & Wegner 1996), quoting broad emission 
lines. 


- VV2000 J093748.7+730158 

It displays a rich absorption system with a complex 
DLA at z = 1.478 and a Lyman limit system near 
z = 2.36 (Rao & Turnshek 2000). 


HS 1002+4400, HS 1049+4033 and HS 1611+4719 
For these three sources, the only references in SIMBAD 
are the report of their discovery in the Hamburg survey 
(Hagen, Engels & Rei mers 1999). Their low-resolution 
spectra (available at http://www.hs.uni-hamburg.de ) 
display rather broad emission lines without very pecu¬ 
liar features. 


- VV2000 J164914.9+530316 


There is not a single explicite reference in SIMBAD ex¬ 
cept the Veron-Cetty & Veron catalogue. As discussed 
in Sect. 3.1, this QSO is weakly radio loud. 


4. Statistical Study 

In order to compare the results of the 1.2 mm studies of 
the z « 2 and z ^ 4 quasars and to search for possible 
correlations between the FIR luminosity (derived from the 
1.2 mm flux density) and the optical luminosity (derived 
from the absolute B-Band magnitude Mb), we made a sta¬ 
tistical study of the present data set and the z ^ 4 QSO 
surveys by Omont et al. (2001) and Carilli et al. (2001a) 
which are based on the PSS and SDSS QSO samples, re¬ 
spectively, and are hereafter referred to as the PSS and 
SDSS samples. The statistical study is based on the sur¬ 
vival analysis method which is the most appropriate to 
handle the censored data (upper or lower limits) in as¬ 
tronomical surveys (Feigelson & Nelson 1985; Isobe et al. 
1986, Isobe & Feigelson 1986). In these methods the values 
of the upper (or lower) limits are formally taken into ac¬ 
count for the statistical analysis. However, the algorithms 
assume that the limiting value is precisely measured, while 
in astronomy the limiting value is rather a specific likeli¬ 
hood upper limit (e.g., Gleser 1992). In this paper, we 
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Fig. 2. The redshift dependence of the 1.2 mm flux den¬ 
sity (S' 250 ) for the 5 w 2 quasars detected in this paper 
and the 2 : > 4 quasars reported by Omont et al. (2001) 
and Carilli et al. (2001a). The open symbols correspond to 
the mean redshift and weighted mean flux density of the 
respective total samples (detections and non-detections) - 
see Table 3. The observational data are compared to the 
flux density expected from a 10^^ Lq source with a spec¬ 
tral energy distribution similar to that of the starburst 
galaxy Arp 220. 


adopt a 3cr (99%) upper limit for the survival analysis 
whereas in the figures of this paper the upper limits are 
shown at the 68% level confidence level (see caption of 

Fig. 0 ). 

4.1. Properties of the z ^ 2 and z > 4 QSOs 

With similar 1.2 mm flux densities, the implied luminosi¬ 
ties, dust masses and star formation rates for the z « 2 
quasars are comparable to the values derived for the SDSS 
and PSS samples. This is due to the fact that at high red- 
shifts for a given FIR spectral energy distribution (SED) 
and luminosity, the (sub)millimeter flux density is nearly 
independent of redshift. At comparable sensitivities, there 
is no obvious difference in the 1.2 mm flux densities of 
the z « 2 and z ^ 4 quasars (note though that the 5250 
flux densities correspond to rest-frame emission at 400 and 
< 240 pm, respectively). This is illustrated in Fig.j^where 
the 1.2 mm flux densities of the z « 2 quasars and of the 
PSS and SDSS quasars are plotted against redshift. 

The statistical properties of the 1.2 mm continuum 
surveys of the present QSO sample (1.8 < z < 2.8) and 
of the PSS and SDSS samples are compared in Table 3. 
Of the 35 z R:! 2 quasars which were observed, 26% were 
detected, which is comparable to the detection fraction of 
29% of the PSS sample (both surveys have comparable 
median r.m.s. values). The 1.2 mm continuum survey of 
the SDSS sample has about twice the depth of the z w 2 
and PSS samples, and a detection fraction of 38%. 


U 

CD 



Fig. 3. The cumulative Mb distribution of the z w 2 QSO 
sample (this paper) compared to the distribution for the 
z > 4 QSO PSS and SDSS samples studied by Omont et 
al. (2001) and Carilli et al. (2001a), respectively. 


A parallel study of z w 2 quasars done at 850 pm by 
Priddey et al. (2002) report a detection fraction of 15% 
for a median r.m.s. error of 2.8 mJy, which is comparable 
to the r.m.s. error of the present study if one adopts a 
typical ratio of 850 pm to 1.2 mm flux densities of 2.5. 

Figure 0 shows the cumulative Mb distributions of the 
observed sources for the z w 2 QSO sample compared to 
those for the PSS and SDSS samples. It illustrates the 
difference in Mb range of the SDSS sample, which has a 
Mb cumulative distribution fainter by about one magni¬ 
tude compared with that of the PSS and z « 2 samples 
(Fig. 0). For the three samples, there is no obvious differ¬ 
ence in the cumulative magnitude distributions between 
the detected sources and the whole sample (not shown in 
Fig.|). 

4.2. The Far-Infrared Luminosity Distributions 

The average 1.2 mm flux density, and hence TfiRj is 
smaller by about a factor two for the SDSS sample com¬ 
pared with the PSS and z w 2 samples (Table 3 and 
Fig. 0). To statistically explore the differences in the FIR 
luminosity distributions between the samples, we applied 
the univariate methods of the survival analysis (e.g., Isobe 
et al. 1986). These statistical tests compute the proba¬ 
bility, p, that two samples derive from the same parent 
distribution. 

Comparing the PSS and SDSS samples, the tests re¬ 
turn probabilities p between 12 and 17%, whereas for the 
PSS and z « 2 samples, p ranges between 25 and 38%, 
and for the SDSS and z « 2 samples, from 1 to 5%. These 
tests indicate that the PSS and z « 2 samples could derive 
from similar parent FIR luminosity distributions, whereas 
the SDSS sample could derive from a different underlying 
distribution. This difference is apparent already from its 
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Table 3. Statistical properties of the 1.2 mm continuum surveys of the present QSO sample (1.8 < z < 2.8) and of 
the z > 4 QSO PSS (Omont et al. 2001) and SDSS (Carilli et al. 2001) samples 


Sample 


Numb. 

Frac. 

(%) 

median r.m.s. 
[mJy] 

(*S’25o) *^5250 

[mJy] 

{Lfir) 
[IQi^ To] 

{Mb) 

z « 2 

Detected 

9 

26 

0.80 

5.24 

0.27 

24.8 

-28.2 


Undetec. 

26 

74 

0.93 

0.71 

0.17 

3.3 

-27.9 


Total 

35 


0.90 

1.97 

0.14 

9.3 

-28.0 

PSS 

Detected 

19 

29 

0.70 

4.52 

0.16 

15.7 

-27.8 


Undetec. 

44 

71 

0.90 

0.40 

0.13 

1.4 

-27.9 


Total 

63 


0.80 

2.07 

0.10 

7.2 

-27.9 

SDSS 

Detected 

15 

38 

0.40 

2.35 

0.11 

8.3 

-27.0 


Undetec. 

25 

62 

0.40 

0.19 

0.08 

0.6 

-27.0 


Total 

40 


0.40 

0.93 

0.06 

3.2 

-27.0 


NOTE - The average flux of the sample is (S' 250 ) = Tj{wiSi)/T,{wi) where the weights are w = l/c^ and 0-3250 ~ [^(^'01”° *■ 
Note that, in the cases of the z « 2 and PSS QSO samples, the average total flux densities and luminosities are biased by the 
lower r.m.s. values of the detections. 

lower average flux densities (Table 3), and is even more 
pronounced when considering the average luminosities. 

That the millimeter fluxes measured for the SDSS 
QSOs are lower than those found for the other samples 
could be related to the lower optical luminosities of the 
SDSS QSO or to some other bias introduced in the source 
selection. Some difference between the three samples could 
also arise from the fact that the PSS and SDSS samples 
result from homogeneous selections, whereas the z « 2 
QSOs are selected from the Veron-Cetty & Veron cata¬ 
logue, which is a compilation of all QSOs discovered by 
various methods. 

The FIR luminosities of the z w 2 QSOs are only 
slightly higher than the luminosities of the z ^ 4 QSOs, 
especially when comparing to the PSS sample which has a 
more similar optical luminosity distribution, a fact which 
indicates that there is no strong FIR luminosity evolution 
for optically luminous QSOs between 2 < z < 5. This re¬ 
sult differs from that obtained by Archibald et al. (2001), 
who reported that the FIR luminosity of radio galaxies is 
a strong function of redshift. 

4.3. Relation between the Far-Infrared and Optical 
Luminosities 

We applied the bivariate methods of the survival analy¬ 
sis to quantify the probability of a correlation between 
the FIR and optical luminosities. For the PSS and SDSS 
samples the tests return high p values indicating that 
there is no correlation between Mb and Tfir- Specifically, 
the Cox regression, the generalized Kendall’s tau and the 
Spearman’s rhcfl tests return values for p of 55, 75 and 
34%, respectively, for the PSS sample, and 77, 86, 94%, 
respectively, for the SDSS sample. For the concatenated 
SDSS and PSS samples p is somewhat lower with 49, 41 
and 15%, respectively. Since p is significantly higher than 

^ The small number of sources, detected and undetected, is 
near the limit of the Spearman’s rho test accuracy 


5%, no correlation between Mb and Tfir can be proved. 
The three bivariate methods applied to the z w 2 QSO 
sample return p = 22, 17 and 10%, respectively, again 
showing no strong correlation between Mb and Lfir. 

The lack of an apparent correlation between the FIR 
and optical luminosities agrees with previous studies 
(Omont et al. 2001; Carilli et al. 2001a; Isaak et al. 2002; 
Bertoldi & Cox 2002; Priddey et al. 2002). However, the 
lower average values of —Mb and S '250 for the the SDSS 
sample relative to the PSS and z ~ 2 samples is sugges¬ 
tive that some correlation could be present, although the 
scatter of the present data is too large to quantify this 
further. QSO samples with a wider range in Mb would be 
required for further studies. 

Some proportionality between the FIR and the optical 
luminosities might well be expected. If the dust is heated 
by the AGN, both the FIR and blue luminosities should 
scale with the bolometric luminosity, and a strong corre¬ 
lation between the two should be seen, and is in fact ob¬ 
served for low-redshift QSOs (Haas et al. 2000). If on the 
other hand young stars are the main heating source of the 
dust, a correlation between the average optical and FIR lu¬ 
minosities could arise from the coeval growth of the central 
black hole and of the stellar spheroid, which is also sug¬ 
gested by the observation of a proportionality of the black 
hole and spheroid masses in local galaxies (e.g., Kormendy 
& Richstone 1995; Magorrian et al. 1998). A correlation 
between the optical and FIR luminosities might be over¬ 
shadowed by a strong scatter due to the uncorrelated tem¬ 
poral variability of star formation and accretion activity. 

5. Discussion 

5.1. Absence of Redshift Evolution of Lpi^ 

The comparison of the 1.2 mm observations of QSOs at 
z ss 2 and z ^ 4 indicates that there is no strong evolu¬ 
tion between 2 < z < 5 of the FIR luminosity function of 
optically bright QSOs. Quasars apparently differ in this re- 
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Fig. 4. Far-infrared luminosities, Lfir, implied by the MAMBO 1.2 mm (250 GHz) flux densities, plotted against the 
redshift z for the 138 QSOs at 1.8 < z < 5.0 reported by Carilli et al. (2001a), Omont et al. (2001) and this paper. 
The sources detected at 1.2 mm are shown as filled symbols, whereas the open symbols represent upper limits to Tfir 
(adopting the same definition as in Fig. for sources not detected at 1.2 mm. The dashed line shows the typical 
r.m.s. flux error of the z « 2 and PSS z > 4 surveys (corresponding to £'250 = 0.8 mJy) and the dotted line shows the 
corresponding 3a detection limit. 


spect from radio galaxies which are FIR-brighter at higher 
redshifts. Archibald et al. (2001) argue that the apparent 
FIR luminosity evolution of radio galaxies they found is 
not a selection effect related to the inevitable youth of the 
known high-redshift radio galaxies. Both their luminos¬ 
ity and morphological evolution lead to the interpretation 
that high-z radio galaxies trace the formation and evolu¬ 
tion of massive elliptical galaxies (van Breugel et al. 1998; 
Lacy et al. 2000; Pentericci et al. 2001). 

It may be premature to claim a true difference between 
the QSOs and high-z radio galaxies. For one, we have yet 
observed only the most optically luminous QSOs, which 
are not representative for the bulk of the QSO emission at 
high redshift, which is due to QSOs with Mb —24, —25. 
Given the locally observed black hole to stellar bulge mass 
ratios, the time-integrated rate of star formation (assum¬ 
ing that is reflected by the FIR emission) vs. black hole 
growth (reflected by the optical emission) must in fact be 
higher by a factor ten compared to what we imply from the 
FIR to optical luminosity ratio of the QSOs in our sam¬ 
ple. This suggests that either the optically less luminous 
QSOs show star formation rates similar to their bright 


peers, or the star formation rate increases in time for the 
entire population - a trend which is not confirmed for our 
QSO sample, and contrary to what is seen for the high-z 
radio galaxies. 

Another important unknown is the effective temper¬ 
ature of the FIR emission - which enters in high order 
when evaluating the FIR luminosity, not even to mention 
a possible evolution of the dust emission properties. 

The lack of an evolution of the apparent FIR bright¬ 
ness (sampled at rest wavelengths 400 /im and 240 /rm at 
redshifts 2 and 4, respectively) of optically bright QSOs 
might well be the cumulative effect of several underlying, 
evolving processes, such as the decline of star formation 
activity along with a drop in effective dust temperature, 
or opposite variations in the contributions of starbust and 
QSO radiation to dust heating. 

Future studies should explore the FIR emission of op¬ 
tically fainter QSOs, and try to measure the temperature 
of the FIR SED to better estimate the FIR luminosity. To 
compare with an independent indicator of star formation 
activity, the molecular line emission could be measured. 
Currently, only the millimeter-brightest QSOs were ob- 
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Fig. 5. Far-infrared luminosities, Afir, implied by the MAMBO 1.2 mm (250 GHz) flux densities as a function of 
rest-frame absolute B-band magnitude, Mb, for the 138 QSOs at 1.8 < z < 5.0 reported by Carilli et al. (2001a), 
Omont et al. (2001) and in this paper (same symbols as in Fig. |4|). 


served, and the few resulting CO line detections (e.g., Cox 
et al. 2002 b) show no redshift dependence of the line to 
continuum flux ratio. 

5.2. Source Counts and Relation between Far-Infrared 
and Optical Luminosities 

From correlated SCUBA and Chandra X-ray observations 
it is estimated that about 10 — 15% of the submillime¬ 
ter background sources contain AGN (Bautz et al. 2000; 
Barger et al. 2001; Blain et al. 2002). Most of these sources 
are not optical quasars, but have optically obscured nuclei. 
Among ^100 blank-field sources detected with SCUBA 
and MAMBO, only two are known to be optical quasars 
(Knudsen et al. 2001; Bertoldi et al. 2000a), 

At millimeter wavelengths the density of background 
sources with S '250 > 2 mJy is observed to be ~ 1000 deg“^ 
(Bertoldi et al. 2000a; Blain et al. 2002). The optically 
very bright quasars we have studied with MAMBO should 
make a negligible contribution to this millimeter back¬ 
ground, as the following simple estimate illustrates. Fan 
et al. ( 2001 a, b, c) infer a number density « 2 deg“^ for 
QSOs with Mb < —26.5 at all redshifts, of which about 
one quarter would be brighter 2 mJy if we assume that 
the millimeter brightness of QSOs is about constant from 


z = 5 — 3.6 (Carilli et al. 2001a; Omont et al. 2001) 
through z w 2 (this work) to the present. Optically very 
bright QSOs targeted in this study should thus only con¬ 
tribute some 0.05% to the millimeter background above 2 
mJy. However, there are ~ 100 times as many QSOs with 
fainter optical luminosities Mb = —23 to —26.5. If these 
had similar millimeter properties as their brighter peers 
(a yet unfounded extrapolation), such QSOs could con¬ 
tribute some 5 % to the millimeter background, a num¬ 
ber marginally consistent with the direct count of such 
objects from the SCUBA and MAMBO surveys and the 
analysis of the Chandra-millimeter associations. However, 
it seems that such counts favour a smaller contribution 
which would imply a slightly smaller QSO FIR luminos¬ 
ity function in the range of Mb« —23 to —25 than their 
brighter counterparts. 

6. Conclusion 

We have presented results of an on-going survey measur¬ 
ing the thermal dust emission of high-redshift quasars. 
After surveying the millimeter continuum of z ^ 4 QSOs 
(Omont et al. 2001; Carilli et al. 2001a), we here re¬ 
port 250 GHz (1.2 mm) observations of optically lu¬ 
minous (Mb < —27.0), radio-quiet quasars at redshift 
1.8 < z < 2 . 8 , corresponding to the peak of quasar ac- 
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tivity. As for the z ^ 4 sources, dust masses of ~ 10® Mq 
and FIR luminosities of Ffir ~ 10^® Lq are implied by 
the observed 1.2 mm flux densities. If the FIR luminos¬ 
ity arises from massive star formation, the inferred star 
formation rate would be w lO®M 0 yr“^. The millimeter 
properties of the 2 « 2 quasars are similar to optically 
comparably bright QSOs at 2 ^ 4. There is no sign of a 
redshift evolution of the FIR luminosity of optically lu¬ 
minous QSOs, in contrast to what was reported for radio 
galaxies by Archibald et al. (2001). 

With a significant sample size in the combined QSO 
samples at redshifts 2 « 2 and 2^4, (43 detections and 95 
non-detections), we find a large scatter and no clear cor¬ 
relation between the millimeter and optical luminosities. 
Lower average values of Mb and S 250 for the SDSS sample 
compared with the PSS sample are suggestive of a gross 
correlation, but the large scatter makes this correlation of 
low statistical significance. 

To better quantify the millimeter properties of high 
redshift QSOs and especially the contribution of star for¬ 
mation to their emission, we need to measure the millime¬ 
ter emission of optically fainter QSOs. We further need to 
constrain the dust temperature through observations at 
submillimeter and far-IR wavelengths, the AGN activity 
through observations at far-IR, radio, and X-ray wave¬ 
lengths, and we must measure the strength of the molecu¬ 
lar and atomic/ionic line emission. Most of these observa¬ 
tions will need to wait for a new generation of sensitive in¬ 
frared/submillimeter observatories: SIRTF, Herschel and 
especially ALMA. 
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